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Results from the ARGO-YBJ Test Experiment 

The ARGO-YBJ Collaboration* 
Presented by M. lacovacci 

'^Dipartimento di Fisica dell'llnivcrsita di Napoli and I.N.F.N., Napoli, Italy 

An RPC carpet covering ~ 10'* (ARGO-YBJ experiment) will be installed in the YangBaJing Laboratory 
(Tibet, P.R. China) at an altitude of 4300 m a.s.L. A test-module of ~ 50 m^ has been put in operation in 
this laboratory and about lO'^ air shower events have been collected. The carpet capability of reconstructing the 
shower features is presented. 



1. INTRODUCTION 

The ARGO-YBJ experiment is under way over 
the next few years at Yangbajing High Alti- 
tude Cosmic Ray Laboratory (4300 m a.s.L, 606 
g/cm^), 90 km North to Lhasa (Tibet, P.R. 
China). The aim of the experiment is the study of 
cosmic rays, mainly cosmic 7-radiation, at an en- 
ergy threshold of ~ 100 GeV , by means of the 
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detection of small size air showers at high al- 
titude. The apparatus consists of a full cover- 
age detector of dimension ~ 71 x 74 realized 
with a single layer of Resistive Plate Counters 
(RPCs). The area surrounding the central detec- 
tor core, up to 100 x lOOm^, consists of a guard 
ring partially (~ 50 %) instrumented with RPCs. 
These outer detector improves the apparatus per- 
formance, enlarging the fiducial area, for the de- 
tection of showers with the core outside the full 
coverage carpet. A lead converter 0.5 cm thick 
will cover uniformly the RPC plane in order to 
increase the number of charged particles by con- 
version of shower photons and to reduce the time 
spread of the shower front. The site location (lon- 
gitude 90° 31' 50" E, latitude 30° 06' 38" N) per- 
mits the monitoring of the Northern hemisphere 
in the dechnation band -10° < d < 70°. 

Such a detector, performing a continuous high 
sensitivity sky survey, complements the narrow 
field of view air Cerenkov telescopes allowing to 
bridge the GeV and TeV energy regions and to 
face a wide range of fundamental issues in Cosmic 
Ray and Astroparticle Physics including 7-ray as- 
tronomy, GRBs physics and the measurement of 
the p/p at TeV energies jl) . Detector assembling 
will start late in 2000 and data taking with the 
first ~ 750 of RPCs in 2001. 

In order to investigate both the RPCs perfor- 
mance at 4300 m a.s.L and the capability of 
the detector to sample the shower front of atmo- 
spheric cascades, during 1998 a full coverage car- 
pet of ^ 50 has been put in operation in the 



2 



YangBaJing (YBJ) Laboratory. Results concern- 
ing both the RPCs behaviour and the features of 
the showers imaged by the carpet are presented. 

2. THE EXPERIMENTAL SET-UP 

The set-up of the detector is an array of 3 x 5 
chambers of area 280 x 112 cmP each, covering 
a total area of ^ 8.6 x 6.0 m^. The active area 
of ~ 45.8 m^, accounting for a dead area due to 
a 7 mm frame closing the chamber edge, corre- 
sponds to a ~ 89% coverage. The RPCs, with a 
2 mm gas gap, are built with bakelite electrode 
plates of volume resistivity in the range (0.5 1) 
■ cm. The RFC signals are picked up by 
means of aluminum strips 3.3 cm wide and 56 
cm long which are glued on a 0.2 mm thick film 
of Poly-Ethylcn-Tcrcftalatc (PET). At the edge 
of the detector the strips are connected to the 
front-end electronics and terminated with 50 fl 
resistors. A grounded aluminum foil is used to 
shield the bottom face of the RPC and an ex- 
tra PET foil, on top of the aluminum, is used as 
a further high voltage insulator. The front-end 
circuit contains 16 discriminators, with about 50 
mV voltage threshold, and provides a FAST-OR 
signal with the same input-to-output delay (10 
ns) for all the channels. This signal is used for 
time measurements and trigger purposes in the 
present test. The 16 strips connected to the same 
front-end board are logically organized in a pad 
of 56 X 56 av? area. Each RFC is therefore sub- 
divided in 10 pads which work like independent 
functional units. The pads are the basic elements 
("pixel") which define the space-time pattern of 
the shower; they give indeed the position and the 
time of each detected hit. At any trigger occur- 
rence the times of all the pads are read-out by 
means of multihit TDCs of 1 ns time bin, op- 
erated in common STOP mode. Since the pad 
signal is shaped to 1.5 fis, only the time profile of 
the earliest particles hitting the pads is imaged. 
The set-up was completed with a small telescope 
consisting of 3 RPCs of 50 x 50 cm^ area with 16 
pick-up strips 3 cm wide connected to front-end 
electronics board similar to the ones used in the 
carpet. The 3 RPCs were overlapped one on the 
other and the triple coincidence of their FAST- 



OR signals was used to define a cosmic ray cross- 
ing the telescope. 

3. RPC PERFORMANCE 

The RPCs were operated in streamer mode as 
foreseen for the final experiment. This mode de- 
livers large amplitude saturated signals and is 
less sensitive than the avalanche or proportional 
mode to electromagnetic noise, to changes in the 
environment conditions and to mechanical de- 
formations of the detector. Three gas mixtures 
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Figure 1. Detection efficiency vs operating volt- 
age for one of the carpet RPCs (•). The same 
curve for a 2 mm gap RPC operating at sea level 
is also reported (o) for comparison. 



were tested which used the same components, 
Argon, Isobutane and TetraFluoroEthane, in dif- 
ferent proportions: TFE/Ar/i-But = 45/45/10; 
60/27/13 and 75/15/10. An higher TFE con- 
centration in place of the Ar concentration in- 
creases the primary ionization thus compensat- 
ing for the 40% reduction caused by the lower 
gas target pressure (600 mbar) and reduces the 



3 



afterpulse probability. The reduction of the Ar- 
gon concentration in favour of TFE results in a 
clear increase of the operating voltage as expected 
from the large quenching action of TFE. Since the 
higher the TFE fraction, the lower is the charge 
delivered in the gas by a single streamer, we de- 
cided to operate the test carpet with the gas mix- 
ture corresponding to the highest TFE fraction, 
in order to extend the dynamic range achievable 
for the analogical read-out. 

Fig. shows the operating efficiency for 4 
ORed pads. The efhciency was measured using 
cosmic ray signals defined by means of an auxil- 
iary telescope placed on top of the carpet. The 
same curve for a 2 mm gap RFC operated at sea 
level is also shown for comparison. The detection 
efficiency vs operating voltage, compared to the 
operation at 606 mbar pressure in YBJ, shows an 
increase of ~ 2.5 kV in operating voltage. The 
effect of small changes in temperature T and pres- 
sure P on the operating voltage can be accounted 
for by rescaling the applied voltage Va according 
to the relationship: 
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where Pq and Tq are arbitrary standard values, 
e.g. 1010 mbar and 293 K respectively for a sea 
level laboratory. This formula predicts, starting 
from the YBJ data, an operating voltage at sea 
level which is considerably smaller than the ex- 
perimental one. However, a good consistency is 
recovered by assuming that, in the ideal gas ap- 
proximation, the parameter which fixes the oper- 
ating voltage is given by gap ■ pressure / temper- 
ature. Measurements on RPCs with different gas 
gap size justify this assumption. 

Fig. 2 also shows that the plateau efficiency 
measured at YBJ is 3 4% lower than at the 
sea level. Although a lower efficiency is expected 
from the smaller number of primary clusters at 
the YBJ pressure, we attribute most of the dif- 
ference to the underestimation of the YBJ effi- 
ciency. At YBJ altitude, indeed, the ratio of the 
cosmic radiation electromagnetic to muon compo- 
nent is ~ 4 times larger than at sea level. A spa- 
tial tracking with redefinition of the carpet track 
downstream would eliminate the contamination 
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Figure 2. Counting rate (full symbols) and oper- 
ating current (open symbols) vs Voltage of one 
RFC of the carpet. Results are presented for 
two different TFE percentages (45%:triangles and 
75%:circles). The rate shown refers to 4 of the 10 
ORed RFC pads. 



from soft particles, giving a more accurate and 
higher efhciency. On the other hand the lower 
efficiency could hardly be explained with the gas 
lower density. In fact, the number of primary 
clusters in the YBJ test, estimated around 9, is 
the same as in the case of some gas, e.g. Ar/i- 
But/CF3Br = 60/37/3, that was frequently used 
at sea level with efficiency of ~ 97 98%. 

The counting rate of 4 ORed pads, together 
with the RFC current, are reported in Fig. ^ vs 
the operating voltage. The results for a percent- 
age of 45% of TFE are also reported for compari- 
son. A rather flat singles counting rate plateau is 
observed at a level of ^ 400 Hz for a single pad 
of area 56 x 56 cm^. 

The time jitter distribution of the pad sig- 
nals was obtained by measuring the delay of the 
FAST-OR signal with respect to RPC2 (the RFC 
in the central position of the trigger telescope) by 
means of a TDC with 1 ns time bin. This distri- 
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Figure 3. Time jitter distribution of 4 pads of 
the carpet. The telescope RFC 2 signal is used as 
common stop. 



bution is shown in Fig. |3| for four pads. The 
average of the standard deviations is 1.42 ns cor- 
responding to a resolution of ~ 1 rts for the single 
RFC if we account for the fact that the distribu- 
tions show the combined jitter of two detectors. 
A more detailed presentation of these results is 
reported in Q. 

4. DATA ANALYSIS OF SHOWER 
EVENTS 

Data were taken either with or without a 0.5 
cm layer of lead on the whole carpet in order to in- 
vestigate the converter effect on multiplicity and 
angular resolution. A trigger based on pad mul- 
tiplicity has been used to collect ~ 10^ shower 
events in April-May 1998. The integral rate as a 
function of the pad multiplicity is shown in Fig. 
^ for showers before and after the lead was in- 
stalled. A comparison at fixed rate indicates an 
increase of pad multiplicity due to the effect of 
the lead of 15 20%, as expected according 
to our simulations. Since in this test the strip 



Figure 4. The integral rate as a function of the 
pad multiplicity. 



read-out is not effective, the actual multiplicity 
is understimed for particle densities higher than 
~ 1 m^^. Data have been corrected to obtain 
the average particle density as a function of the 
recorded multiplicity (Fig. |[). 

4.1. Event reconstruction 

Due to the reduced size of the detector, the 
time profile of the shower front sampled by the 
carpet is expected to exhibit a planar shape. In 
this approximation the expected particle arrival 
time is a linear function of the position. We per- 
forme an optimized reconstruction procedure of 
the shower direction as follows: 

• Unweighted plane fit to hits for each event 
by minimization of the function 

= \ y^^{lxi + myi + nZi-c(ti-tQ)y{2) 

i 

The sum includes all counters with a time 
signal ti, c is the light velocity, (xi,yi,Zi) 
are the ith center pad position coordinates. 
The parameters of the fit are the time offset 
<o and the I, m direction cosines. 
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Figure 5. The particle density as a function of 
the measured multiplicity. 



• Rejection of outlying points by means of a 
K ■ a cut and iteration of the fit until all 
points verify this condition. If the remain- 
ing hits number is < 5 the event is rejected. 
Here a is the standard deviation of the dis- 
tribution. 

This procedure is rather fast because it makes 
use only of analytic formulae. No a priori infor- 
mation about shower features is required. The 
actual value of K could depend on the features 
of the reconstructed showers as well as on the ex- 
perimental conditions (pad dimension, shaping of 
the signals, multihit capability, etc.). By choos- 
ing K = 2.5 about 8% of the hits are rejected. 
Increasing K should increase the number of dis- 
carded hits without a significant improvement of 
the x^- The distribution exhibits a long tail 
more pronounced for low multiplicity events. In 
these high events the time hits are consider- 
ably spread, by far more than expected from the 
detector time resolution. They can be attributed 
to the sampling of a portion of the shower affected 
by large time fluctuations. Events with reduced 
> 30 ns^ (~ 10% of the total) have been dis- 



carded and not used in the following analysis. 

4.2. The zenith angle distribution 

Since the shower size is not measured, events 
have been classified according to the observed 
particle density. The zenith angle distribution 
of the integral density spectrum is expected to 
follow an exponential behaviour 

Ii>P,0)=I{>p,0)-e-^^'''''-^^ (3) 

where xq is the vertical depth. The attenua- 
tion mean free path Aatt of EAS is related to 
the absorption length Ae of EAS particles by 
Aatt = Ae/7 where 7 is the index of the integral 
density spectrum. The angular distribution of 
events with density greater than 3.0 particles/m^ 
(corresponding to Nhit{0) = 76 • cos6) is shown in 
Fig. ^ . The data can be fit out to ~ 55° with 
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Figure 6. The zenith angle distribution. 



an exp[—a ■ {secO — 1)] law. The parameter a is 
found to be 4.88 ± 0.45, so that A^tt = (124 ± 11) 
g/cm?, in excellent agreement with previous re- 
sults concerning the density spectrum of shower 
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particles ||3|]. For angles greater than 55° a de- 
viation from this law is observed. Misrecon- 
structed events, showers locally produced in the 
walls of the sourrounding buildings and horizontal 
air showers could contribute to these large angle 
events. This result is substantially unchanged for 
densities ranging from 1 to 5 particles/m^. 

Data confirm that the shape of the angular dis- 
tribution is dominated by the physical effect of 
atmospheric absorption and that no relevant in- 
strumental effects or selection biases have been 
introduced. 

4.3. The density spectrum 

The differential density spectrum measured in 
three different intervals of zenith angles is shown 
in Fig. . The angular bin size is A0 = 5°, 
so that the difference of atmospheric depth inside 
the source bin is < 40 5 /cm^ (about one radiation 
length) for zenith angles up to 35°. A smaller bin 
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Figure 7. Differential density spectrum for three 
intervals of zenith angle. 



size has been used to analyze data up to 55°. The 
spectra refer to the depth corresponding to the 



mean zenith angle of the showers recorded inside 
each interval. The shapes of the spectra are very 
similar in the density range 2-^8 particles/m^ not 
affected by threshold or saturation effects. 

The index 7^ of the density spectrum as a func- 
tion of the zenith angle is shown in Fig. ^. A 
weighted mean up to 9 < 45° gives 2.33 ± 0.03. 
This value is fully consistent with the ones ob- 
tained in experiments at high altitude at com- 
parable particle densities MonteCarlo simu- 
lations show that present data concern energies 
ranging from about 1 TeV (quasi-vertical show- 
ers) to a few hundred TeV {0 ^ 45°). Accord- 
ingly, the measured jd can be considered a rea- 
sonable estimate of the slope of the density spec- 
trum of showers initiated by primaries (mainly 
protons) of energies ~ 10^^-;-10^^ eV. From these 
data the absorption length of EAS particles Ag 
is found to be 165 ± 18 g/cw? consistent with 
results of similar measurements at small shower 
sizes HJ. Data recorded at large zenith angles 
(> 45°) are due to showers developed through at- 
mospheric depths greater than 860 gjcm?. Thus, 
a substantial contribution is expected from show- 
ers originated by primaries of energies beyond the 
knee of the cosmic ray spectrum. This could ex- 
plain the increasing of the slope 7(j with zenith 
angle, as shown in Fig. ||. 

4.4. Shovi^er front thickness 

The time distribution of the shower hits with 
respect to the fitted plane is shown in Fig. |9| for 
two different multiplicity ranges. Quasi-vertical 
showers {Q < 15°) have been selected. The dis- 
tribution of time residuals exhibits, as expected, 
a long tail due to time fluctuations and to the 
curved profile of the shower front, more pro- 
nounced for low multiplicity. The width of these 
distributions is related to the time thickness of 
the shower front. Since the position of the shower 
core is not reconstructed, the experimental result 
concerns a time thickness averaged on different 
radial distances. 

For low multiplicity events with a mean number 
of particles per pad < 1, the spread a represents a 
reasonable measurement of the shower thickness. 
For high multiplicity events {Nhu > 100), where 
the mean number of particles hitting one pad is 
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Figure 8. The index 7^ of the density spectrum 
as a function of the zenith angle. 



Figure 9. Distribution of time residuals for 
events with different pad multiplicity (all chan- 
nel added). 



> 1, this a does not express the thickness of the 
shower disc without any bias. In fact, the time 
residual distribution is related to the fluctuations 
of the first particle detected by each pad. Taking 
into account the total detector resolution of 1.3 
ns (RFC intrinsic jitter, strip length, electron- 
ics time resolution) the time jitter of the earliest 
particles in high multiplicity events (> 100 hits) 
is estimated 1 ns. The averaged shower front 
thickness of low multiplicity events (particle den- 
sity < 1 rn~^) is found ^ 4.4 ns. 

4.5. Angular resolution 

The angular resolution of the carpet has been 
estimated by dividing the detector into two inde- 
pendent sub-arrays (" odd pads" and " even pads" ) 
and comparing the two reconstructed shower di- 
rections. These two sub-arrays overlap spatially 
so that they are affected by the same shower cur- 
vature. Events with total number of hits Nhu 
have been selected according to the constraint 
Nodd ^ Neven ^ Nhit/2. The distribution of the 
even-odd angle difference A6'eo is shown in Fig. 
fol for events in different multiplicity ranges and 



9 < 55°. These distributions narrow, as expected, 
with increasing shower size. 

The effect of the lead sheet on the angular res- 
olution can be appreciated in Fig. |l^, where the 
median M^g^^ of the distribution of A9eo as a 
function of pad multiplicity, for showers recon- 
structed before and after the lead was added, is 
shown. The improvement of the angular resolu- 
tion is a factor 1.4 for Nhu = 50 and decreases 
with increasing multiplicity. Assuming that the 
angular resolution function for the entire array is 
Gaussian, its standard deviation cr is given by |^] 

(Tg ~ i iY7x2 ■ '^^^^ chessboard method yields a 
standard deviation ag ^ 2.1° for events with a 
pad multiplicity > 100. 

5. CONCLUSIONS 

A carpet of ~ 50 of RPCs has been put in 
operation at the Yangbajing Laboratory in order 
to study the high altitude performance of RPCs 
and the detector capability of imaging with high 
granularity a small portion of the EAS disc, in 
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Figure 10. Even-odd angle difference distribution Figure 11. Median of A6eo distribution as a func- 
for events with different pad multiplicity. tion of pad multiplicity. 



view of an enlarged use in Tibet (ARGO-YBJ 



experiment). ^^^^ .^^^2) 350 

The results of this test confirm that RPCs can ^ ' 

be operated efBciently (> 95%) to sample air 
showers at high altitude with excellent time res- 
olution 1 ns). The analysis of data collected 
with a shower trigger suggest that the RPCs car- 
pet capability of reconstructing the shower fea- 
tures is consistent with expectations. As a con- 
clusion, the overall results of the test look well 
promising for future operation of the full detec- 
tor. 



Ill (1956). 

Alexandrcas D.E. et al., Nucl. Instr. Meth. 
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